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We present the first evidence of the decay B — > r v T , using 414 fb _1 of data collected at the T(4S) resonance 
with the Belle detector at the KEKB asymmetric-energy e + e _ collider. Events are tagged by fully reconstructing 
one of the B mesons in hadronic modes. We detect the signal with a significance of 4.0 standard deviations includ- 
ing systematics, and measure the branching fraction to be B{B~ — > t~u t ) = (l-06^Q'2g (stat)^'^ (syst)) X 10 -4 . 
We also report results based on 1.86fb _1 data collected by the Belle detector at the T(5S) resonance. Several 
exclusive B a decays B s — > D B *^ ir~ (p~) and B s — » J/ijitj>(rf) are studied. The B a meson production is found to 
proceed predominantly through the creation of B*B* pairs. Upper limits on B s — > K + K~ , B a — » (jyy, B s —* 77 



and B s 



D&+D. 



(*)- 



decays are also reported. 



1. Introduction 

In the Standard Model (SM), the purely leptonic 
decay B~ — > t~v t proceeds via annihilation of b and 
u quarks to a W~ boson. It provides a direct determi- 
nation of the product of the B meson decay constant 
Jb and the magnitude of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element \V u b\- The branching 
fraction is given by 



GpmBrn^. 



57T \ m h 

where Gf is the Fermi coupling constant, rag and 
m T are the B and r masses, respectively, and tb is 
the B~ lifetime [1]. Physics beyond the SM, such as 
supersymmetry or two-Higgs doublet models, could 
modify B(B~ — > t~D t ) through the introduction of 
a charged Higgs boson [2]. Purely leptonic B decays 
have not been observed in past experiments. The most 
stringent upper limit on B~ — > t~D t comes from the 
BaBar experiment: B(B~ -> t~z> t ) < 2.6xl0" 4 (90% 
C.L.) [3]. 

The possibility to study decays of B s at very high 
luminosity e + e~ colliders running at the energy of the 
T(5S f ) resonance has been discussed in several the- 
oretical papers [4, 5]. The first data at the T(55) 
were taken many years ago at CESR [6-8], but the 
collected data sample was not enough to observe a 
B s signal. In 2003, the CLEO experiment collected 
0.42 ftT 1 at the T(55) and observed some evidence 
for B s meson production in both inclusive and exclu- 
sive modes. However, simple calculations assuming 
an approximate SU(3) symmetry indicate that many 
interesting B s measurements require a data sample of 
at least 20 pb _1 , which can be collected by B Facto- 
ries in the future. To test the experimental feasibility 
of such measurements, a data sample of 1.86 fb _1 was 
recently taken with the Belle detector at the center- 
of-mass (CM) energy corresponding to the mass of the 
T(5S f ) resonance. This data sample is more than four 
times larger than the CLEO dataset at the T(55). 

The Belle detector is a large-solid-anglc magnetic 



spectrometer consisting of a silicon vertex detector, 
a 50-layer central drift chamber (CDC), a system of 
aerogel threshold Cerenkov counters (ACC), time-of- 
flight scintillation counters (TOF) , and an electromag- 
netic calorimeter comprised of CsI(Tl) crystals (ECL) 
located inside a superconducting solenoid coil that 
provides a 1.5 T magnetic field. An iron flux-return 
located outside of the coil is instrumented to identify 
if£ and muons. The detector is described in detail 
elsewhere [9]. 



2. Evidence of the Purely Leptonic 
Decay B~ — > t~v> t 

We use a 414 fb _1 data sample containing 447 x 10 6 
B meson pairs collected with the Belle detector at 
the KEKB asymmetric-energy e + e~ (3.5 on 8 GeV) 
collider operating at the T(45) resonance (-y/s = 10.58 
GeV). 

We use a detailed Monte Carlo (MC) simulation, 
which fully describes the detector geometry and re- 
sponse based on GEANT [10], to determine the signal 
selection efficiency and to study the background. In 
order to reproduce effects of beam background, data 
taken with random triggers for each run period are 
overlaid on simulated events. The B~ — > t~D t signal 
decay is generated by the EvtGen package [11]. To 
model the background from e + e~ — > BB and contin- 
uum qq (q — u, d, s, c) production processes, large BB 
and qq MC samples corresponding to about twice the 
data sample are used. We also use MC samples for 
rare B decay processes, such as charmless hadronic, 
radiative, electroweak decays and b — > u semileptonic 
decays. 

We fully reconstruct one of the B mesons in the 
event, referred to hereafter as the tag side (B tag ), and 
compare properties of the remaining particle(s), re- 
ferred to as the signal side (B s i s ), to those expected 
for signal and background. The method allows us to 
suppress strongly the combinatorial background from 
both BB and continuum events. In order to avoid ex- 
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perimental bias, the signal region in data is not looked 
at until the event selection criteria are finalized. 

The Stag candidates are reconstructed in the fol- 
lowing decay modes: B + — > XJ(*)°7r + , 
D^°a+ andjDM°E>i* )+ . The D mesons are recon- 
structed as D° -> K+tt-, K+ 7T 7T°, K + TT TT + 7T , 
K^n°, K%tt-tt+, K^tt-jt+tt and and the 

Df mesons are reconstructed as Df — > i£gif + and 
K + K~ir + . The E>*° and D* + mesons are recon- 
structed in D*° -> D 7T°,!D 7, and £>*+ -> E>+ 7 
modes. The selection of -B tag candidates is based 

on the beam-constrained mass Mb c = V ^bcam — Pb 
and the energy difference AE = Eb — Et, cam . Here, 
Eb and are the reconstructed energy and momen- 
tum of the i?tag candidate in the e + e~ center-of-mass 
system, and Eb oa m is the beam energy in the CM 
frame. The selection criteria for Stag are defined as 
M bc > 5.27 GeV/c 2 and -80 MeV < AE < 60 MeV. 
If an event has multiple B tag candidates, we choose the 
one having the smallest x 2 based on deviations from 
the nominal values of AE, the D candidate mass, and 
the D* — D mass difference if applicable. 

In the events where a -Btag is reconstructed, we 
search for decays of B slg into a r and a neutrino. 
Candidate events are required to have one or three 
charged track(s) on the signal side with the total 
charge being opposite to that of Btag- The r lep- 
ton is identified in the five decay modes, ^jl~v^i> t , 
e~v e v T , n~v T , Tt~'K v T and ir~ ir + ir~ v T , which taken 
together correspond to 81% of all r decays [1]. The 
muon, electron and charged pion candidates are se- 
lected based on information from particle identifica- 
tion devices. The leptons are selected with require- 
ments that have efficiencies greater than 90% for both 
muons and electrons in the momentum region above 
1.2 GeV/c, and misidentification rates of less than 
0.2%(1.5%) for electrons (muons) in the same mo- 
mentum region. Kaon candidates are rejected for all 
charged tracks on the signal side. The tt° candidates 
are reconstructed by requiring the invariant mass of 
two 7's to satisfy |M 77 - m w o\ < 20 MeV/c 2 . For all 
modes except t~ — ► ir^ir ^, we reject events with ir° 
mesons on the signal side. All the selection criteria 
have been optimized to achieve the highest sensitivity 
in MC. 

The most powerful variable for separating signal 
and background is the remaining energy in the ECL, 
denoted as Eecl, which is sum of the energy of pho- 
tons that are not associated with either the Stag or 
the 7T° candidate from the r~ — > -K^n a v T decay. For 
signal events, Eecl must be either zero or small value 
arising from beam background hits, therefore, signal 
events peak at low Eecl- On the other hand, back- 
ground events are distributed toward higher Eecl due 
to the contribution from additional neutral clusters. 

The -Eecl signal region is optimized for each t de- 
cay mode based on the MC simulation, and is de- 
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Figure 1: -Eecl distribution for the control sample of 
doubly tagged events, where one B is fully reconstructed 
and the other B is reconstructed as B~ —* D*°l~v. The 
dots with errors indicate the data. The solid histogram 
represents the distribution as deduced background from 
BB MC (B + B~ + B°B°), and the dashed histogram 
shows the contribution from B°B° events. 



fined by -Eecl < 0.2 GeV for the ii~v^v T , e~v e v T and 
•n~v T modes, and -Eecl < 0.3 GeV for the tt~7t°v t and 
7r~7r + 7r~^ r modes. The -Eecl sideband region is de- 
fined by 0.4 GeV < -Eecl < I- 2 GeV for the [Tv^Vt, 
e~v e v T and tt~v t modes, and by 0.45 GeV < Eecl < 
1.2 GeV for the 7r _ 7r°i/ r and 7r~7r + 7r~j/ T modes. Ta- 
ble I shows the number of events found in the sideband 
region for data (-V°j^|) and for the background MC 
simulation (iV s ^) scaled to the equivalent integrated 
luminosity in data. Their good agreement for each r 
decay mode indicates the validity of the background 
MC simulation. Table I also shows the number of the 
background MC events in the signal region (N^p). 

In order to validate the Eecl simulation, we use a 
control sample of events (double tagged events) , where 
the -Bt a g is fully reconstructed as described above 
and -B s jg is reconstructed in the decay chain, B~ — ► 
D* £~D (D*° -> D°tt°), followed by D° -> K~7T + 
or K~tt~tt + tt + where £ is a muon or electron. The 
sources affecting the Eecl distribution in the control 
sample are similar to those affecting the Eecl distri- 
bution in the signal MC simulation. Figure 1 shows 
the Eecl distribution in the control sample for data 
and the MC simulation scaled to equivalent integrated 
luminosity in data. Their agreement demonstrates the 
validity of the Eecl simulation in the signal MC. 

After finalizing the signal selection criteria, the sig- 
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Figure 2: -Eecl distributions in the data after all 
selection requirements except the one on -Eecl have been 
applied. The data and background MC samples are 
represented by the points with error bars and the solid 
histogram, respectively. The solid curve shows the result 
of the fit with the sum of the signal shape (dashed) and 
background shape (dotted). 



nal region is examined. Figure 2 shows the obtained 
-Eecl distribution when all r decay modes are com- 
bined. One can see a significant excess of events in the 
Eecl signal region below Eecl < 0.25 GeV. Table I 
shows the number of events observed in the signal re- 
gion (N },s) for each r decay mode. For the events in 
the signal region, we verify that the distributions of 
the event selection variables other than Eecl, such as 
Mbc and p m i ss , are consistent with the sum of the sig- 
nal and background distributions expected from MC. 

We deduce the final results by fitting the obtained 
Eecl distributions to the sum of the expected sig- 
nal and background shapes. Probability density func- 
tions (PDFs) for the signal / s (Eecl) and for the back- 
ground ^(Eecl) are constructed for each r decay 
mode from the MC simulation. The signal PDF is 
modeled as the sum of a Gaussian function, centered 
at Eecl = 0, and an exponential function. The back- 
ground PDF, as determined from the MC simulation, 
is parametrized by a second-order polynomial. The 
PDFs are combined into an extended likelihood func- 
tion, 



C 



-(n s +n b ) 



where Ej is the Eecl in the ith event, N is the total 
number of events in the data, and n s and are the 
signal yield and background yield to be determined 
by the fit. To combine likelihood functions of the five 
decay modes, we multiply the likelihood functions to 
produce the combined likelihood (£ com = 11^=1 
The results are listed in Table I. The number of signal 
events in the signal region deduced from the fit (N s ) 
is 2I.2I5 7 when all t decay modes are combined. Ta- 
ble I also gives the number of background events in 
the signal region deduced from the fit (N^), which is 
consistent with the expectation from the background 
MC simulation (iV s ¥ c ). 

The branching fractions are calculated as B = 
N s /(2 ■ e ■ Nb+b-) where Nb+b- is the number of 
T(4S) — ► B + B~ events, assumed to be half of the 
number of produced B meson pairs. The efficiency is 
defined as e = e tag x e scl , where e tag is the tag re- 
construction efficiency for events with B~ — > t~v t 
decays on the signal side, determined by MC to be 
0.136 ± 0.001(stat)%, and e scl is the event selection 
efficiency listed in Table I, as determined by the ratio 
of the number of events surviving all of the selection 
criteria including the r decay branching fractions over 
the number of fully reconstructed E ± . When all r de- 
cay modes are combined we obtain a branching frac- 
tion of (1-06^0 28) x 10~ 4 . The branching fraction for 
each r decay mode is consistent within error as shown 
in Table I. 

Systematic errors for the measured branching frac- 
tion are associated with the uncertainties in the num- 
ber of B + B~ , signal yields and efficiencies. The total 
fractional uncertainty of the combined measurement 



is 

be 



+20.5 

-24. 



q%, and we measure the branching fraction to 



B(B~ 



v T ) = (1.06l°;^( S tat)+£i(syst)) x 10 



+0.22/ 



nw s ( £ i)+^( £ i)). (2) 



The significance is 4. Oct when all r decay modes are 
combined, where the significance is defined as £ = 
y/—2 ln(£o/£max), where £ m ax and Co denote the 
maximum likelihood value and likelihood value ob- 
tained assuming zero signal events, respectively. Here 
the likelihood function from the fit is convolved with a 
Gaussian systematic error function in order to include 
the systematic uncertainty in the signal yield. 



3. Results from the T(55) Engineering 
Run 

We use a data sample of 1.86 fb _1 taken at the 
T(5S I ) energy of ~ 10869 MeV. The experimental con- 
ditions of data taking at T(55) are identical to that 
for T(45) or continuum running. The data sample of 
3.67 fb _1 taken in the continuum at an energy of 60 
MeV below the T(45) was also used in this analysis 
for comparison. 



fpcp06_114 



4 



Flavor Physics and CP Violation Conference, Vancouver, 2006 









N MC 


N ohs 


N s 




e scl (%) 


B(10- 4 ) 


£ 




96 


94.2 ± 8.0 


9.4 ±2.6 


13 


5.4±|J 


g.i±g-; 


8.88 ± 0.05 


i-oi±g-S 


2. Oct 




93 


89.6 ± 8.0 


8.6 ±2.3 


12 


3.9+H 


9-2 + °o1 


8.18 ±0.05 


0.79 +0 £ 


1.3(7 


■k~v t 


43 


41.3 ±6.2 


4.7 ± 1.7 


9 


3.4t;-« 


4.ol° ;? 


5.79 ± 0.04 


0.96lS:S 


1.9(7 


■K~ -K°V T 


21 


23.3 ±4.7 


5.9 ±1.9 


11 




4.2l° ; 3 3 


8.32 ± 0.08 


i 9 o+0.77 
1 -^°-0.B3 


2.3cr 


1T~ 7T + 7T~ h> T 


21 


18.5 ±4.1 


4.2 ± 1.6 


9 


o i +3.1 
°- 1 -2.6 


n 7 +0.3 
°- ' -0.2 


1.75 ± 0.03 


9 QQ+ 3 01 
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274 


266.9 ± 14.3 


32.8 ±4.6 


54 


21.21?:? 


30.2l°j 


32.92 ±0.12 


1 06+ ' 34 
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Table I The number of observed events in data in the sideband region (N°^ B C ), number of background MC events in the 
sideband region (iV^) and the signal region (A^ c ), number of observed events in data in the signal region (iVobs), 
number of signal (N s ) and background (Nh) in the signal region determined by the fit, signal selection efficiencies (e scl ), 
extracted branching fraction (£>) for B~ — ► t~v t . The listed errors are statistical only. The last column gives the 
significance of the signal including the systematic uncertainty in the signal yield (E). 



The B s mesons are produced at the T(5S) through 
the intermediate B S B S , B*B S , B S B* or B*B* pair 
production channels, where B* decays to B s j. These 
intermediate channels can be distinguished kinemati- 
cally and their production ratios can be obtained from 
the reconstruction of exclusive B s decays. To improve 
the statistical significance of our exclusive B s signal, 
we combined the six modes B s — > £>+7r~, B s — > 
D*+7r-, B s -» D+p~, B s -» D*+p-, B s -» J/ij,<f> and 
B s — » J/V>?7 , which have large reconstruction efficien- 
cies and are described by unsuppressed conventional 
tree diagrams. 

Six conventional B s decays to -D+7r~, D+p~, 
D* + 7r - , D* + p~, J/ip<j) and J/tpr] final states and four 
rare B s decays to K + K~, (jry, 77 and D s D s 
final states are reconstructed. The signals can be 
observed using two variables: the energy difference 
AE = E™ — -Ebcam an d beam-constrained mass 

M bc = ^(£ b C I) 2 - (P™) 2 ; E b M and p™ are the 
energy and momentum of the B s candidate in the 
CM system and E™ m is the CM beam energy. The 
B s mesons can be produced at the T(55) energy via 
the intermediate e + e~ — ► Bg*^ channels, with 
B* — > B s j. The B s signal regions in Mb c and AE 
are separated for different intermediate channels. 

After all selections, the dominant background is 
from e + e~ — ► qq continuum events (q — u, d, s, or c). 
The distribution of data in Mb c and AE for the B s — > 
Dfjr~ decay mode is shown in Figure 3a. Three Z)+ 
decay modes, <f)ir + , -ft"* -/^ and iYgif+j are used to 
reconstruct B s candidates. Nine events are observed 
within the B s signal ellipsoidal region corresponding 
to B*B* pair production channel. Only one event is 
observed in the B s signal region for B*B S + B S B* 
channels, and no events are observed for B S B S chan- 
nel. Background outside the signal regions is small 
and corresponds to 0.1 event for any of three signal 
regions. The inclusive studies at the T(55) found that 

92, 000 ± 7, 900 ± 23, 500 B^B^ pairs are contained 
within that 1.86 fb _1 T(5S) data sample. Using 
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Figure 3: The M bc and AE scatter plot for 
B s -> Dfir-(a), B s -> D* s + iv-(b) and B s -> D ( s * )+ p~ (c) 
decay modes are shown, where Df meson is 
reconstructed in the Df — ► <j)ir + , Df — ► K*°K + and 
Df — > KgK + decay modes. Also Mbc and AE scatter 
plot (d) is shown for the B s — ► J/ip<f> and B s — » J/tpV- 



this value, we measure the branching fraction to be 
B(B 8 -> D+7T-) = (0.65 ± 0.21 ± 0.19)%. 

The Mbc and A-E scatter plots are also obtained 
for the B s D* s +tt- (Figure 3b) and B s D { * )+ p~ 
(Figure 3c) decay modes. Again, three Df decay 
modes, <f)ir + , K*°K + and KgK + , are used to recon- 
struct B s candidates. The numbers of events within 
the B s signal region for the B* S B* pair production 
channel are 4 for B s — > D* + tt~ decay and 7 for 

£> s — ► Z?s decay. 

The scatter plot in Mb c and AE for the B s — > J/ -00 
and -B s — > J/ipn decays is shown in Figure 3d. One of 
the observed B s — > J/V"^ candidates is reconstructed 
in the J/^> — > A* + / i_ mode and one in the J/^ — > e+e - 
mode. These two candidates correspond roughly to a 
branching fraction of <~ 1 x 10 -3 , in agreement with 
expectations. 

The B s and B* masses can be extracted from the 
Mbc fits in the B*B* channel. The Mb c distribution 
for this channel (Figure 4a) is obtained choosing can- 
didates within the -0.08 < AE < -0.02 MeV range. 
The distribution, shown in Figure 4a, is fitted by the 
sum of a Gaussian to describe the signal and the AR- 
GUS function to describe the background. The fit 
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Figure 4: The B* (a) and _B S (b) mass distributions for 
events within the -0.08 < AE < -0.02 MeV interval, 
corresponding to the B*B* channel. 



yields the mass value M(B*) = 5418±1 MeV/c 2 . The 
observed width of the B s signal is 3.6 ± 0.6 MeV/c 2 
and agrees with the value obtained from the MC sim- 
ulation. Using events from the B*B* channel we can 
obtain also the B s mass (Figure 4b). The distribu- 
tion shown in Figure 4b is fitted to the sum of a 
Gaussian and the ARGUS function. The fit yields the 
B s mass M(B S ) = 5370 ± 1 ± 3 MeV/c 2 and width 
a{B s ) = 3.6 ±0.6 MeV/c 2 . 

Additionally, we searched for several B s rare decays 
for the first time: the penguin decay B s — > K + K~ , 
the electromagnetic penguin decay B s — > (fry, and the 
intrinsic penguin decay B s — > 77. We also searched 
for the tree decay B s — ► Dl*^ + , which is not 
yet observed and is of special interest because the 
states are expected to be dominantly CP 
eigenstates. Although the branching fractions for 



these decays are expected to be too low for observation 
in this analysis, we obtained upper limits (Table II). 



Decay mode 


Yield 
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Eff. 


upper limit 




events 


events 
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(xHT 4 ) 


B s -> K+K~ 


2 


0.14 


9.5 


3.4 


B s 07 


1 


0.15 


5.9 


4.1 


B s -> 77 





0.5 


20.0 


0.56 


B s -> D+D- 





0.02 


0.020 


710 


Bs D*+D- 


1 


0.01 


0.0090 


1270 


B s D*s+Dr 





< 0.01 


0.0052 


2730 



Table II The number of events in the signal region, 
estimated background events and upper limits on the 
90% confidence level for B 3 -> K + K~ , B s -> 07, 
Bs 77 and B s — ► D s *^ + D s *^ decay modes. 
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